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Для дифференциальных уравнений, описывающих одномерное стационарное течение с переходом через скорость звука, пред- 
лагаются зависимости в правых частях этих уравнений, позволяющие получить аналитические выражения для параметров газа 
в канале постоянного сечения. 


1. Введение 


При решении стационарных газодинамических 
задач возникает проблема апробации разработанных 
программ путем сравнения с точными решениями. 
Для нестационарных течений таких решений доста- 
точно много [1], в то время как для стационарных - 
значительно меньше. Прежде всего, это уравнения, 
описывающие течение от источника (стока), и транс- 
цендентное уравнение, описывающее одномерное 
распределение числа Маха вдоль канала переменного 
сечения при изоэнтальпическом, изоэнтропическом 
течении идеального совершенного газа [2]: 


тіп(^4) 


= М 


у + 1 


2 + (у-1)М- 


(Г+1)/(2(Г-1» 


( 1 ) 


где: М=іНр/(уР) - число Маха; р, II, Р, А, у -плот- 
ность, скорость, давление, площадь поперечного 
сечения сопла, показатель адиабаты газа. Большин- 
ство других решений базируется на (1). Так, в [3] на- 
хождение параметров двухфазного течения сводит- 
ся к решению уравнения (1) путем использования 
гипотезы о законе отставания частиц и введением 
некоторого эффективного показателя адиабаты. 

В данной работе предлагаются зависимости в 
правых частях уравнения движения и энергии, при 
которых дифференциальные уравнения имеют 
точное решение в виде трансцендентного уравне- 
ния для распределения числа Маха вдоль канала 
постоянного сечения, и явные выражения для ос- 


тальных параметров газа. Более простые соотноше- 
ния приведены в [4] . 


2. Исходные уравнения 

Рассмотрим одномерные стационарные уравне- 
ния для идеального совершенного газа в канале 
постоянного сечения: 

с/рГ/_ 0 сІ(рЦ 2 +Р) _ с сІрЦН _ {2 

ах ах ах 


где: // ' , г ■ Щ - полная энтальпия; х - про- 

у-Ір 2 

дольная координата, принадлежащая рассматрива- 
емой области [х„;х к ] ; С, и С 2 - правые части уравне- 
ний движения и энергии. 

Полагается, что: у=сошІ; заданы граничные ус- 
ловия на входе в сопло в виде Н=Н(х а ), 5=5(х а ), 
С 1 (х а )=С 2 (х а )= 0; внутри рассматриваемой области 
число Маха монотонно возрастает от дозвуковой до 
сверхзвуковой величины, и существует только одна 
точка х,, в которой М= 1 . 

Переходя к переменным Н, 5, М, N и используя 
вместо уравнения неразрывности его интеграл, пе- 
репишем уравнения (2) в следующем виде: 

(іБ = [ с і г 1 (7-0 (3) 
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где: 


7Ѵ = 
7+1 


(5) 


с__Р о _ Г+1 о _І 

р 1 ” ° 2 ( 7 - 1 )’ 7’ / 1-г' 


Старые переменные выражаются через новые с 
использованием формул: 
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, Р = 8р г . (6) 


С учетов вышеупомянутых ограничений на чис- 
ло Маха можно показать [5], что для N получается 
необходимое и достаточное условие существования 
в точке х, единственного минимума: 


тіп(ТѴ) = 7Ѵ(х,), 
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Произведем обезразмеривание отнесением х - к 
ширине канала, II - к критической скорости II,, р- к 
критической плотности р., Р - к произведению р,К 2 . 
Тогда Н(х а )=Н 0 и 8(х а )=8 0 . Подставляя последние два 
выражения в (5) получим граничное условие для N 
Щх а ) = 1 . ( 8 ) 


3. Точные решения 

Пусть заданы х,, Ъ ь Ь ь причем х а <х.<х ь , 0 <Ь 2 , 
0<Ь 2 <\. Рассматривая простейшие зависимости, 
представим 5 в виде линейной функции, удовлет- 
воряющей граничным условиям и условию не убы- 
вания энтропии, а N - в виде квадратичной функ- 
ции, удовлетворяющей (7, 8): 

Б = Б 0 (Ь І (х-х,) + \), (9) 
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Тогда число Маха находится из трансцендент- 
ного уравнения (4), точное решение для //получа- 
ется из (5): 


Н = Н 0 



а остальные параметры находятся из (6). Диффе- 
ренцируя (5) по "х" и, используя (3), получим соот- 
ношения для С,, С 2 , которые можно использовать в 
(2) для численного решения: 
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причем /V, 8, /V', 8' в (11) берутся из (9), (10). 
Штрих обозначает производную по "х. Очевидно, 
что если в (9) положить /,=0, то этому будет соотве- 
тствовать случай 5=5' 0 =соп8і. 

Для Я=// =сопм и ( 10) число Маха находится из 
(4), точное решение для 8 получается из (5): 

8 = Б 0 мК 


остальные параметры из (6). Естественно, что в 
этом случае энтропийная функция на определен- 
ном участке может убывать. Правые части для ур. 
(2) получаются из ур. (3) и (5): 

С 2 =0, С=Р*. (12) 

где N и Я' берутся из (10). 

Очевидно, что в предложенных соотношениях 
распределение числа Маха, определяемое из (4), 
зависит только от величин х а , х,, Ь 2 в (10). 

4. Численные расчеты 

На рис. 1 приведено распределение числа Маха 
вдоль канала, полученное из ур. (4) и ( 10) при зна- 
чениях 7=1,4, х=-4, х.=0, х ь =2, />2={0,1; 0,5; 0,9}. 

Для предложенных правых частей ур. (11, 12) 
был проведен ряд расчетов методом установления 
по явной схеме Мак-Кормака [6] для нестационар- 
ного аналога ур. (2): 

ф + фП =0 дрЦ | д(рЦ 2 + Р) _ с 
ді дх ’ ді дх 1 ’ 

д(рн-Р) | дрцн _ с 
ді дх 2 ’ 

число точек сетки при этом равнялось 40. Началь- 
ное распределение, показанное на рис. 1, находи- 
лось из соотношений: 

Ч = 0,18(х + 2) +1, Р = Б оР \ 


что соответствует начальному положению точки 
х,=-2. 



Рис. 1. Распределение числа Маха вдоль канала 


В процессе установления при использовании 
правых частей ур. (11, 12) и некоторого значения Ъ 2 , 
распределение числа Маха вдоль канала будет схо- 
диться к соответствующей кривой на рис. 1, при 
этом положение точки х, будет сходиться к 0. 

На рис. 2 и рис. 3 показано положение самой 
левой из точек х, в процессе установления. 
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Технические науки 



О 400 число итераций §00 

Рис. 2. Положение х> в процессе установления при Ьу=0,4 



Рис. 3. Положение х, в процессе установления при Н=сопзІ 

Видно, что при Ъ\= 0,4 и правых частях (11) амп- 
литуда колебаний больше, чем при правых частях 
( 12). Как видно из рисунков, во всех расчетах поло- 
жение точки х, сходится к точному значению, рав- 


ному 0, причем во всех случаях остается только од- 
на точка х,. 

Несмотря на то, что метод установления демон- 
стрирует сходимость к точному решению и предель- 
ное распределение числа Маха вдоль канала, полу- 
ченное при 2000 итераций (рис. 4) отличается от 
точного, изображенного на рис. 1 в третьем знаке 
после запятой, на промежуточных итерациях вполне 
возможно нарушение монотонности и появление 
нескольких точек х,, что демонстрирует рис. 4. На 
рис. 2 и рис. 3 резкие скачки при й,=0,9 как раз свя- 
заны с наличием нескольких точек х, и сопровожда- 
ющейся при этом перестройкой течения. 



Рис. 4. Распределение числа Маха вдоль канала для разных 
итераций. I — номер итерации 


5. Заключение 

Предложены аналитические формулы, позво- 
ляющие находить распределение параметров газа с 
переходом через скорость звука в канале постоян- 
ного сечения. Эти формулы целесообразно исполь- 
зовать при тестировании соответствующих числен- 
ных методов. 
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